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ABSTRACT 

Two scenarios have been proposed for evolution of star forming cores: grav- 
itational fragmentation of larger structures and coalescence of smaller entities 
which are formed from some instabilities. Here, we turn our attention to the lat- 
ter idea to investigate the evolution of observed low-mass condensations (LMCs) 
in the cores of molecular clouds. For this purpose, we implement the evolution of 
the observed LMCs of Taurus molecular cloud 1 (TMC-1). The core is modeled 
as a contracting cylinder with randomly spawned condensations in the middle 
region around its axis. For advancing bodies in their trajectories, we represent 
the acceleration of a particular LMC in terms of a fourth-order polynomial using 
the predictor-corrector scheme. Whenever two LMCs collide, they are assumed 
to be merged in one large condensation containing all the masses of the two pro- 
genitors. Implementations of many computer experiments with a wide variety 
of the free parameters show that the LMCs merge to form star-forming regions 
in the core. The results show that the total mechanical energy of the core in- 
creases by time, and its rate of increasing decreases by facilitating the merger. 
Finally, the mass spectrum index and goodness-of-fit are determined with 50% 
error in the number of mass points. The results show that the goodness-of-fit 
will be refined at the end of simulations, and the mass spectrum index inclines to 
the observed values for the moderate mass objects. The simulations show that 
the TMC-1 turns about 40% of its mass into cluster of dynamically unstable 
protostellar cores. In general, we suggest that the future of LMCs in a core of 
molecular cloud is merger to convert about half of its initial masses into a cluster 
of gravitationally unstable protostellar cores. 

Subject headings: ISM: clouds - ISM: evolution - stars: formation - Methods: 
numerical. 
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1. Introduction 

Star formation appears to occur exclusively within the molecular phases of the inter- 
stellar medium (see, e.g., Klessen et al. 2009). Molecular clouds exhibit extremely complex 
hierarchical structure, with column densities that vary by many orders of magnitude. The 
smallest divisions are nominated as protostellar cores, with masses of a few solar mass or less 
and sizes of few tenth parsecs or less (e.g., Stahler and Palla 2004). The advent of sensitive 
tools and advanced observational methods in the last decade has made it possible to perform 
high-resolution surveys to interior of the individual cores (e.g., Tafalla 2008). In this way, 
many of embedded condensed objects within each core have been revealed that we call them 
low-mass condensations (LMCs). Although a lot of these have recently found by the Spitzer 
c2d project (e.g., Lee et al. 2009), but we document this research on the pioneer report of 
Peng et al. (1998, hereafter PLVKL) who presented observational data on the small-scale 
structure of the Taurus Molecular Cloud 1 (TMC-1) in the regime of 0.02 — 0.04pc and 
0.04 — O.O6M . This preference is because of many small groups and associations of star for- 
mation in the Taurus molecular cloud that provides it as an important test-place for theories 
and models (Hartman 2002). 

The formation mechanisms of LMCs are not well understood. In fact, two schemes 
have been proposed: gravitational fragmentation of larger structures and merger of smaller 
fragments. Preserving the latter mechanism requires some suitable theoretical models for 
the formation of the smaller entities in the cloud. One model that may be regarded is the 
effect of turbulence, although estimations of the level of turbulence in the molecular cores 
give subsonic values (e.g., Myers 1998, Lee et al. 2004, Andre et al. 2007). As another 
idea, van Loo et al. (2007) implemented a two-dimensional simulation with the idea that 
the mechanism causing substructures is the same as the mechanism causing core formation 
(i.e., hydromagnetic waves). The appeared substructures in their simulation are accurately 
conformed to the observed substructures of core D of TMC-1. Although, the dynamical 
models like turbulence or hydromagnetic waves can be considered as sources of substructures 
in dense cores, but we must not neglect some non-dynamical processes and instabilities such 
as thermal instability (e.g., Fukue and Kamaya 2007, Nejad-Asghar 2007, Nejad-Asghar and 
Molteni 2008, Banerjee et al. 2009, Nejad-Asghar and Soltani 2009). In any case, formation 
of a system of LMCs in each protostellar core is a controversial issue. In this research, we 
have nothing to do with the formation issue of LMCs, rather we attempt to model their 
future. Since, our aim is to present a realistic model for future evolution of LMCs, we use 
the data of PLVKL for LMCs in the core D of TMC-1. The results may approximately be 
applicable to all cores with a system of LMCs. 

The high-resolution investigation in the southeastern part of the TMC-1, which was 
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done by PLVKL, covers an elongated 8' x 8' area in the sky, or 0.32pc at the distance of 
140pc. The emission in this region originates in three narrow components centered on Local 
Standard of Rest (LSR) velocities of ~ 5.7, ~ 5.9, and ~ 6.1km.s _1 . These components 
each represent a separate cylindrical feature elongated along the ridge of TMC-1. Among 
the three velocity components, PLVKL identified a total of 45 LMCs with masses in the 
range of 0.04 — O.O6M and sizes of 0.02 — 0.04pc. PLVKL suggested that thermal and 
turbulent pressures of the ambient gas of the core may contribute to keep these LMCs 
pressure bond so that the dispersion time of LMCs is comparable to the molecular cloud's 
lifetime ~ 10 6 yr (e.g., Tan et al. 2006). Typical inflow velocities in the protostellar cores are 
of order 0.05 — O.lkm.s -1 (Lee et al. 1999), which is comparable to the results of PLVKL 
for the random motion of order 0.05km.s _1 among the LMCs. We suggest the gravitational 
force between LMCs can influence the future of them. On the other words, over the lifetime 
of LMCs, they are more likely to be subject to merge because of subsonic collisions with 
each other rather than to travel intact along the core. Investigating this future is the task 
of this research. 

In fact, modeling future evolution of LMCs in protostellar cores adequately requires 
the accurate and simultaneous treatment of many different physical processes. Here, we try 
to overcome the computing limitations by using an adapted N-body method, in which the 
magnetic field, shock waves, heating and cooling processes are ignored. The effect of ambient 
gas on motion of each LMC is assumed to be only a drag force as outlined by Ostriker (1999). 
Whenever two LMCs approach each other, the result of interaction is determined according 
to prescribed rules, which have been numerically examined by many authors (e.g., Lattanzio 
and Henriksen 1988, Gittins et al. 2003, Kitsionas and Whitworth 2007). The simulations, in 
general, show that in subsonic collisions with small impact parameter, the outcome is always 
a merger. Burkert and Alves (2009) have recently used this idea to simulate the future of 
the two LMCs in the protostellar core Barnard 68. Their simulation shows that merger of 
these two sub-cores will lead to gravitational instability. Thus, merger of LMCs might in 
general play an important role in triggering star formation and shaping the molecular core 
mass distribution. 

It seems that using a suitable gravitational N-body method for evolution of LMCs in the 
core D of TMC-1, and using a parametric method for considering the merger of them, can 
give us some acceptable progress for future of the observed LMCs in the molecular cores. For 
this purpose, setup of computer experiments with initial positions and velocities of LMCs in 
the core D of TMC-1 are introduced in section 2. Advancing the LMCs in their trajectories, 
considering of merger, and the results of simulations are presented in section 3. Finally, 
section 4 devotes to a summary and conclusion. 
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2. Setup of computer experiments 

We choose the units of length and time as U[ = 3 x 10 14 m = O.Olpc and u t = 1.5 x 10 4 yr, 
respectively, therefore the unit of velocity is u v = 0.67km.s _1 . The mass unit is chosen as 
u m = 2 x 10 30 kg = 1M Q so that the gravitational constant is G = 1. Each LMC in the core D 
of TMC-1 is assumed to be a sphere with average radius R = 3.39 x lQ~ 2 {b ma jb m i n ) l l 2 ui where 
a distance of 140pc is used, and b ma j and 6 m j„ are major and minor axes in arc-seconds, 
respectively, which are listed in Table 1 of PLVKL. The radii of the 45 LMCs and their 
masses at the present time, t = 0, are shown in Fig. 1. If we assume that the LMCs are 
pressure confined, self-gravitating isothermal gas spheres in hydrostatic equilibrium, we can 
approximately adopt the relation R 18.7M10~ M between the radius and mass of them in 
units of U[ and u m , respectively (see Appendix A). 

The region that analyzed by PLVKL is an 8' x 8' area of the TMC-1 which is equivalent 
to 0.32pc for the cloud at distance of 140pc. The distribution of integrated intensity of the 
three components b (blue), m (middle), and r (red) are mapped in Figure 2 of PLVKL. 
Overall, maps show an elongated appearance, which may represent a cylindrical structure 
for the core. In this way, we consider a parent core in a cylindrical geometry with length 
Ahui and diameter 15^. The LMCs of blue (b) and middle (m) components are distributed 
along the length of cylinder, while the LMCs of red (r) component are approximately placed 
in three quarter of the length. We divide the realm of each LMC by some subdivisions along 
the length of the parent core, and then its center is randomly placed into these subdivisions. 
The azimuth angle for position of each LMC is randomly determined between to 2tt, and its 
axial distance is randomly chosen between fi n Ro and f ou tRo where R = 7.5?/; is the radius of 
cylindrical core. According to the best matches to the depicted maps of Figure 2 of PLVKL, 
the fractions fi n and f out are obtained approximately equal to 0.4 and 0.7, respectively (see 
Fig. 2). 

As shown in Table 1 and 2 of PLVKL, the speed of LMCs have a small dispersion 
about the mean with a standard deviation of 0.03km.s _1 = 0.045u v . On the other hand, 
Lee et al. (1999) show that the typical inflow velocities in the molecular cores are of order 
0.05 — 0.1km. s _1 . Thus, here we assume the speed of each LMC, v, to be randomly generated 
in the range of 0.03 — O.lkm.s -1 . The velocity vector has two components: a radial part v r 
toward the center of the core, and a transverse component v±. As depicted in Fig. 3, the 
velocity in the rc'yV-coordinate is given by 

v' = i'v± cos <p' + j'v± sin ip' + k'v r , (1) 

where tp' is the azimuth angle. 

We introduce two free parameters for velocity vector: (1) the transverse velocity param- 
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eter a so that v± = av and v r = \j\ — a 2 v, and (2) the inclination parameter f3 so that the 
azimuth angle <p' is randomly produced between ±/3tt/2. In this way, the velocity vector v' 
for each LMC is determined whenever the values of these two free parameters are appointed. 
Since, the coordinate x'y'z' is obtained via two successive rotations of xyz, 

( — cos #0 sin 9 \ / cos <p sin ip \ 
1 ) ( -sinp cosp ) , (2) 
-sinfl -cosfl / V 1 / 

the velocity vector v in rcyz-coordinate can be determined by 

v = [R y ,{9 + Tr)R z {<p)]- 1 v'. (3) 



Each LMC within the aggregate of bodies in the parent core experiences an acceleration 
that arises from two sources: (1) the gravitational attractions of other LMCs, 



a (0 

grav 



and (2) a drag force from ambient gas, 



4t 9 = J i G2P i M ^ ( 5 ) 



s 



where p and c s are the mass density and sound speed of the ambient gas in the parent core, 
respectively (Ostriker 1999). 

The description of the problem is now completed by specified accelerations for all the 
LMCs. Thus, the position of the bodies can be advanced by an accurate N-body method. 
Whenever the impact parameter of two colliding LMCs, with radii Ri and Rj, is less than 
(Ri + Rj)/^Y, they are assumed to be merged in a one larger condensation with radius ~ 
18.7M1CT M (see Appendix A), containing all the mass of the two progenitors (M = Mj+M,). 
The collision is assumed to be perfectly inelastic so that the position and velocity after merge 
are given by 

(6) 
(7) 

Mi + Mj ' K ' 

respectively. Note that the properties of the models in each computer experiment are com- 
pletely controlled by three key parameters a, /3, and 7. 



and 



MjVj + MjTj 
Mi + Mj 



Mm + MjVj 
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3. Results of simulations 

Determining the new positions of the LMCs at a somewhat advanced time is deceptively 
simple in appearance, but a slight carelessness can lead to a complete change of long-term be- 
havior. Here, we use the method of Aarseth (2003) for advancing LMCs in their trajectories. 
This method is a version of a predictor-corrector scheme, which pivots on representing the 
acceleration of a particular body of index i at time t in terms of a fourth-order polynomial 
based on knowledge of the acceleration at four previous times in the past, to, ti, t2, and t^, 
with t being the most recent, 

a(t) = (((D 4 (t - t 3 ) + D 3 )(t - t 2 ) + D 2 )(t - h) + Di)(t - t ) + a , (8) 

where a = a. grav + a.d rag is the acceleration at time t - Using compact notation, the first 
three divided differences are defined by 

L>k[to,tk\ = , {k = 1, 2, 3), (9) 

£o — tfc 

where Do = a and square brackets refer to the appropriate time intervals. The term D4 is 
defined similarly as 

d 4 = v&M-vM (10) 

t — £3 

which the predictor corrector method is used to obtain quantities at the advanced time t 
that information is not initially available. 

When the acceleration of each LMC is calculated from (8), one can immediately integrate 
with respect to time obtaining its velocity and then again to obtain the position. The key 
idea underlying to do is determining a suitable time-step. A simple method for estimating 
the time-step for each LMC involves the distance and the relative velocity with respect to 
its nearest neighbor. Here, we adopt the more sophisticated relation 

/ |a||a( 2 )| + |aW| 2 \ 1/2 
At i = V\ rjrnrj^i , rj^T, , (H) 



where 77 is a dimensionless accuracy parameter (rs 0.01 in our computer experiments), and 
the superscript in parentheses is the order of the derivative (i.e., a^) = da/dt, and so on). 

The properties of models in the computer experiments are completely controlled by 
three key parameters: transverse velocity component < a < 1, axial inclination of the 
velocity vector < (3 < 2, and merger parameter 7^4. Simulations with different values of 
these free parameters show that in any case, the LMCs merge to form more massive larger 
bodies. The results, as expected, show that decreasing of the three parameters a, /3, and 7 
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leads to facilitating the merger of LMCs and formation of star-forming regions in the parent 
core. With typical values of a = 0.1, ft = 1.0, and 7 = 4.0, the number of LMCs, N, and 
the mean mass, J^Mj/AT, at different simulation time are shown in Fig. 4. 

In a perfectly inelastic collision between two LMCs which leads to merger, the kinetic 
energy dissipates and the two-body potential energy approaches to zero. Thus, in each tow- 
body merger, the total mechanical energy increases. This fact is shown in Fig. 5 for three 
values of the merger parameter 7, in which the parameters a and ft are chosen equal to 0.1 
and 1.0, respectively. As can be seen in this figure, the total mechanical energy of a system of 
LMCs, which are initially gravitationally bond, increases by time until it approaches to zero 
at a critical time t c . Afterwards, the system will be gravitationally unbound. Fortunately, 
this critical time is less than the lifetime of LMCs so that the results of our computer 
experiments, in this simplest implementation, are reliably fruitful at t < t c . The rate of 
increasing of the total mechanical energy of a system of LMCs, in consequence of merger, 
decreases by decreasing of the parameters a, j3, and 7. On the other words, the critical 
time t c increases by facilitating the merger of LMCs in the computer experiments. The 
computer experiments show that the effect of the parameters a and ft on the critical time 
t c is negligible, although the results indicate that decreasing of the transverse parameter a 
and inclination parameter j3 can facilitate the merger of LMCs and leads to little increasing 
of critical time t c . 

The mass spectrum of clumps in giant molecular clouds appear to be well described by 
a power law, 

£«M->, (12) 

with the index being in the range —1.8 < a — 1 < —1.3 (e.g., Kramer et al. 1998). When 
considering the protostellar cores, the mass function is well described by a double power law 
fit, following approximately a — 1 r; —2.5 above ~ O.5M and a — 1 r; —1.5 below (Motte 
et al. 1998, Testi and Sargent 1998, Johnstone et al. 2001, Johnstone et al. 2006, Alves et 
al. 2007). There can also be found a power law mass function in the low-regime of brown 
dwarfs (Thies and Kroupa 2008). Overall, the canonical mass function with confidence 
found so far is a power law function with universality hypothesis that it constitutes the 
parent distribution of all stellar populations (Kroupa 2008). 

In the low-mass regime of LMCs, PLVKL excluded the lowest and highest mass points 
to study the mass distribution of these small entities by using a power law fit, logA^ = 
alogM + const., with a m —0.45. They indicated that the corresponding mass spectrum 
index a — 1 = —1.45 is in reasonable agreement with those obtained by observational data 
analysis of stellar populations. Here, we assume that there is approximately 50% error in the 
number of LMCs (A(logiV) ~ | log(1.5)|), and the mass spectrum index a — 1 is evaluated in 
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the simulation time. The results of this linear fit and its goodness-of-fit (Press et al. 1992) 



where Q is the incomplete gamma function and x 2 is chi-square, are given in Fig. 6. This 
figure shows that the goodness-of-fit will be refined at the end of simulation, and the mass 
spectrum index inclines to the expected values for the moderate mass populations. Fig. 7 
shows the histogram of the mass distribution of LMCs at present time, and at the simulation 
time t ~ 1.5 x 10 5 yr for a = 0.1, = 1.0, and 7 = 4. At present time, there is only one 
unstable collapsing LMC, while at the end of simulation, there are five unstable LMCs. 
Thus, the simulation shows that the core D of TMC-1 will be a cluster of five gravitationally 
unstable protostellar cores with masses ~3x 0.8 + 2 x 0.7 = 3.8M . 



Star formation is a great theme in astronomy today, and special attention is attributed 
to the formation of low-mass stellar groups. Large ground-based telescopes and currently 
active satellite observatories have revealed many low-mass condensations (LMCs) in the 
molecular cores. It seems that the collisions and merger of these LMCs may lead to formation 
of moderate mass objects and eventually to formation of small groups of stars within the 
parent core. Here, we made a try to investigate the future of these LMCs. For this purpose, 
as an excellent test-place, we applied the observed substructure fragments in the core D of 
TMC-1, which their radii and masses at the present time are shown in Fig. 1. 

Initial positions of LMCs in the cylindrical parent core are depicted in Fig. 2, in which 
their initial axial distances are randomly chosen between 0.4i? and 0.7-R where R is the 
radius of cylinder. The observations show that the LMCs have infall velocities in the sub- 
sonic regime. Here, we assumed that there is two free parameters for the velocity vector of 
LMCs: transverse component parameter a and axial inclination parameter /3. The radial and 
transverse components of the initial velocity vector of each LMC are schematically shown in 
Fig. 3. The acceleration (gravitational attraction and drag force) of LMCs are represented 
in terms of the fourth-order polynomial, and we used a merger parameter 7 to state the 
coalescence of them. 

There are three key parameters a, /3, and 7 which control the properties of a model in 
the computer experiments. The results of computer experiments with a wide variety of these 
parameters show that overall the LMCs merger to form more massive larger bodies. As a 
direct result of merge, the number of LMCs decreases by time, and their mean mass increases 
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4. Summary and conclusion 
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as shown in Fig. 4 for a typical model. As another result of merger, the total mechanical 
energy of a system of LMCs increases by time until it approach to zero at a critical time t c . 
This critical time increases by facilitating the merger of LMCs which arose by decreasing of 
the parameters a, /3, and 7 (Fig. 5). 

At last, the mass distribution of LMCs are evaluated in the simulation time, and the 
results of a power law fit and its goodness-of-fit, with 50% error in the number of mass points, 
are shown in Fig. 6. On the whole, the results show that a system of LMCs in a low-mass 
protostellar core merge to form small groups of star-forming regions. The simulations show 
that the core D of TMC-1 turns about 40% of its mass into cluster of dynamically unstable 
protostellar cores. In general, we suggest that the future of LMCs in a core of molecular cloud 
is merger to covert about half of its initial masses into a cluster of gravitationally unstable 
protostellar cores. Clearly, over the collapse of each unstable protostellar core in the cluster, 
the mass conversion drops to few percent, consistent with the overall star formation efficiency 
of molecular clouds in the Galaxy. Although, the results of this research are deduced from 
simulation of LMCs in the core D of the TMC-1, but can reliably be applicable for future of 
all molecular cores with a system of LMCs. 
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A. Radius-mass relation for LMCs 

We suppose that LMCs are pressure confined, self-gravitating isothermal gas spheres in 
hydrostatic equilibrium, which conform the Lane-Emden equation, 

1 d ( 2 dif)\ 

where the dependent variable tp = 4>/c 2 s is defined as the ratio of gravitational potential to 
the isothermal sound speed, and £ = (4:irGp c / '(? s ) l l 2 r is the nondimensional length where p c is 
the central density at r = 0. The equation (Al) can be solved numerically by the boundary 
conditions ^(0) = ^'(O) = 0, then the density profile can be obtained by p = p c exp(—tp). 
Since we assume that each LMC is confined by an external pressure p ext , its edge with radius 
imax can be determined by the condition p(£ max ) = p ext /c 2 s . In this way, the nondimensional 



(Al) 
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mass m = (plxtG^/c^M of each LMC is given by 

The Fig. 8 shows the density contrast p c j p(£ ma x) and radius of pressure-bounded isother- 
mal spheres as a function of nondimensional mass. Masses smaller than the critical value 
m = 1.18 are stable while ones with ^ max > 6.5 (or equivalently p c j p(£ ma x) > 14.1) are 
gravitationally unstable so that may collapse. Here, we approximately use the linear fit- 
ted relations to express the numerical results of stable states of Fig. 8 (in the range of 
< m < 1.18) as follows: 

imax = 5.51m, log 10 ( Pc ] = 0.97m. (A3) 

\PKKmax); 

There is a variety range of densities reported in the TMC-1. For instance, around 
the cyanopolyyne peak, Schloerb et al. (1983) and Bujarrabal et al. (1981) found molecular 
hydrogen densities of 5 — 10 x 10 4 and 1 — 3x 10 4 cm~ 3 , respectively, both from multi-transition 
studies of HC3N. Pratap et al. (1997) found that over most of the TMC-1 ridge the density 
was relatively uniform with a value of approximately 6 x 10 4 cm~ 3 . PLVKL used the CCS 
column densities to distinguish LMCs with an H 2 densities of ~ a few xl0 4 cm~ 3 . Here, we 
apply an appropriate mean value of ~ 5 x 10 4 cm~ 3 for the ambient gas of the TMC-1, thus 
we have 

and 

M = 0.5 ( ^ ( - ni trh ) V2 m (A5) 



^ioKy V 5xl ° 4cm ~ 3 , 

for radius and mass of each LMC in units of U[ and u m , respectively. Eliminating the m 
between (A4) and (A5) by using the linear approximations (A3), and choosing T = 10K and 
n(£,max) = 5 x 10 4 cm~ 3 for TMC-1, we approximately have 

R « 18.7M10~ M (A6) 

In this way, the dimensional radius can directly be obtained via equation (A6) for merged 
stable LMCs with total mass Mi + Mj < 0.5 x 1.18. Thus, not only the dimensional radius 
of the merged LMCs are obtained, but also its condition for collapse can be checked. 
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Fig. 1. — The radii of the 45 LMCs in the core D of TMC-1 and their masses at the present 
time t = are shown as the black squares. The curve is the function R = 18.7M10 _M which 
is explained in the Appendix A. 
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Fig. 2. — Initial positions of LMCs in the cylindrical core D of TMC-1. They are grouped 
by their LSR velocity components, with designations by b (blue), m (middle), and r (red) 
referring to the 5.7, 5.9, and 6.1km.s _1 , respectively. The cylindrical parent core is rotated 
in this figure so that it can be easily compared to the Figure 2 of PLVKL. Here, the fractions 
fin and f out are assumed to be 0.4 and 0.7, respectively. 
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Fig. 3. — Schematic diagram of the radial and transverse components of the initial velocity 
vector of each LMC. The origin of x'yV-coordinate is placed on the LMC, and the origin of 
rryz-coordinate is assumed to be at the center of the core. 
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Fig. 4. — (a) Number of LMCs, and (b) the mean mass of them versus the simulation time, 
for a = 0.1, 13 = 1.0, and 7 = 4. 
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Fig. 5. — Total mechanical energy of LMCs versus simulation time for three typical values 
of merger parameter 7. The parameters a and (3 are 0.1 and 1.0, respectively. Clearly, 
decreasing the parameter 7 facilitates the merger of LMCs. 
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Fig. 6. — (a) Mass spectrum index versus simulation time with 7 = 4, a = 0.1 and ft = 1.0. 
(b) The goodness-of-fit with assumption of 50% error in the number of LMCs. 
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Fig. 7. — Histogram of the mass distribution of LMCs at present time (top panel), and at 
the simulation time t ~ 1.5 x 10 5 yr (bottom panel) for a = 0.1, = 1.0, and 7 = 4. 
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Fig. 8. — The density contrast p c j ' p(£ m ax) and the nondimensional radius of pressure-bounded 
isothermal spheres as a function of nondimensional mass. The dashed lines are approximately 
linear relations to the numerical results for the gravitationally stable cases. 



